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ACCURACY OF ESTIMATING THE I4AS3ES OF PHOBOS AND 


DEIMOS FROM MULTIPLE VIKING 0R3ITFR ENCOUNTERS 


Robert H. Tolson and Mary L. Mason 


This paper addresses the problem of estimating 
the masses of Phobos and Peimos from Doppler 
and onboard optical measurements during the 
Viking extended mission. A Kalman filter ic 
used to analyze the effects of gravitational 
uncertainties and non-gravitaticnal accelera- 
tions. These accelerations destroy the dynami- 
cal integrity of the orbit and multi-batch or 
limited memory filtering i3 preferred to single 
batch processing. Optical tracking is essential 
to improve the relative orbit geometry. The 
masses can be determined to about 10£ and 255? 
respectively for Phobos and Deimoe, assuming 
satellite densities of about 3 gr/cm^. 


INTRODUCTION 

Two Viking spacecraft will be placed into orbits about Mars in 1 1 .« siimmer of 
1976. Each spacecraft, will separate into un orbitcr (VO) and a lander (VL). 
Upon separation, the VL will perform a propulsive maneuver, causing the 
spacecraft to enter the atmosphere and soft land on the surface of Kars. 

The orbiters will perform experiments to support the lunders and act as 
relay stations to return VL scientific data to the Earth. 


i 
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The spacecraft orbital parameters are largely determined by the landing sites 
on Mars. For the nominal landing si i the orbital geometry is shown in 
Fig. 1. Both orbiters have period vnchronous with the rotational period of 
Mars and periapsis altitudes of 150o km. The orientation elements are given 
in July for the A spacecraft and August for the B spacecraft. The 
precession of the VO-A spacecraft due to Murs planetary oblateness provides a 

H • 1500 km VO-A: I • 34° 

S - SUBSOtAR POINT H P *J?700 km VOB: I • 70° 

A -LANDING Sin a 


B • LAND INC S 111 PERIOD • 24.61 hr 



Figure 1. Inertial orbits of the Viking orbiters. 

unique opportunity to explore the natural satellites Phobos and Deimos from 
distances measured in tens of kilometers 1 . The geometric aspects of this 
opportunity are shown in Fig. 2. The plant* of the figure is the Mars 
equatorial plane. The points where the spacecraft ascends and descends 
through the equatorial plane are indicated for a number of orbits. Initially, 
the descending intersection point is between the orbital paths of the natural 
satellites and the ascending intersection point is inside the Phobos orbit. 
However, nodal regression and apsidal precession cause both of these points 
to move toward the Phobos orbital path. In January and March the spacecraft 
will pass vithin 10 km of the Phobos orbit. As precession continues, the 
ascending intersection point moves outward so that in late 1977 VO-A will 
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pass within 50 km of the Deimos orbital path. Thus, with proper phasing, the 
spacecraft can be maneuvered to pass very close to both natural satellites. 

Since phasing can be performed over relatively long time periods, it can be 
accomplished for very small propulsive coot. For the January and March 



Figure 2. Fhobos and Del nos inertial encounter 
geometry for VO-A 


encounters, it is also possible to establish u 1 to i ratio between the 
spacecraft and Phobos orbital perio : for about PC m/s of propulsive expen- 
diture^, which is well within the spacecraft capability. Thus, repeated 
encounters at close distances are possible. The Deimos encounter opportunity 
requires more propulsive expenditure to maximize the number - of close encounters, 
because the Deimos orbital period (30.3 hrs) i3 substantially dif'*» - <nt from 
the nominal VO orbital period. However, other commensurable ptricu rutius, 
such as 5 to 1*, can be established for a few meters per second. 



The final encounter geometry will be designed to accomplish a number of 
scientific objectives, among which is the determination of the natural 
satellite masses. The mass of Phobos is a significant geophysical parameter, 
since it can be related directly to internal properties of Mars. In addition, 
if the volumes of the natural satellites can be determined, knowledge cf the 
masses can be used to infer density ani hence, composition. It is expected 

that a combination of Mariner 9 and Viking imagery dita will produce satellite 
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volumes with an accuracy of about 10* . Thus, mass determination to less than 
lOt is desirable. 

This paper addresses the problem of determining the mas.ies of the natural 
satellites from the VO-A tracking data und from imagery data. This parametric 
analysis of the masses of Phobos and Deimos considers various approach 
geometries, dynamic noise levels, and data processing methods. 

ERROR ANALYSIS 

Dynan c Model . The dynamic model for an error analysis does not demand the 
fidelity of a data processing model. To analvze close encounters of Phobos or 
Deimos it is only necessary that the spacecraft and satellite orbital motions 
include the dominant perturbative effects. Perturbations to the two body 
orbital motion are due to the .oon-central part of the Martian gravitational 
field, third body forces, radiation pressure, drag, ana g.s ...eaks. The 
dominant gravitational perturbation is the oblateneso ,.e J 0 term, which 
to first order produces the apsidal precession and non 'egression discus: ?d 
above. Tesseral harmonics, third body forces, etc., produce variations in the 
orbit which will not substantially effect the relative geometry. Thus, the 
spacecraft orbital motion is modeled as a uniformly precessing conic according 
to the classical formulas: 
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The orbital elements, the secular rates, and ether physical parameters for 

Phobos and Deimos used in this analysis are given in Table 1. The orbital 
3 U 

parameters and mean diameters wore determined frem Mariner 9 optical data. 

3 5 6 

More detailed orbital representations are available ’ ; however, a first 

order model is more than adequate for this analysis. The gravi tat tonal 

3 

constant assunes a density of about 3 gr/cm' . 

Observational Model . The Viking spacecraft has two video imaging systems 
which have a small amount of image overlap and parallel optical axes. Each 
image is composed of an array of about 1000 * 1000 picture elements or pixels. 
These systems are calibrated with star images during interplanetary cruise 
to measure image distortion and system alignment. This calibration assures 
an accuracy of all subsequent pictures or picture pairs to the pixel level of 
about 0.029 milliradians . 

The optimum utilization of this system for determin'^g the orbits of the 
natural satellites is to simultaneously image the satellite with one system 

Table 1 


PHYSICAL PARAMETERS OF PHOBOS AND DEIMOS ON NOVEMBER 11 , 1971 



Phobos 

Deimos 

a, am 

9378.53+. 01 

2314 58.911. 03 

e 

0.0150+.0001 

0.0008+. 0001 

Mo, deg 

311.8 ?+. 36 

232.6+7. 

i, deg 

1.0U+.01 

2.79+.02 

0 ), deg 

?69 • 9+ • 9 

235.6+7. 

«, dog 

100. 5+. 8 

10.91.2 

n, deg/day 

11 28 . It 069 

285. 1*4 38 

Ui 

.87148 

.0361 

• 

n 

-.U 37 I 4 

-.0181 

• • 

n + w + deg/day 

11?8. 81*1*3+. 0001 

28 5 . 16181. 0001 

Mea 1 Diameter, km 

23.0 

12.8 

Gravitational Constant 
km 3 /sec 2 

2 x 10" 3 

2 x IQ" 1 * 
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and the adjacent star background with the other . The center of figure of the 
natural catelllte relative to the star background is then determined by fitting 
an ellipsoidal shape to the optical image. The center of the resulting 
ellipsoid is then assumed to be the center of mass. The two major sources 
of error are the basic oixel resolution and the center of macs determination 
within the image. To account for density inhomogenity and surface irregular- 
ities, it is assumed that this latter process has ar. error of ', 0 % of the 
diameter of the natural satellite. Thus the total optical data error is 
modeled as the HSS of the 0.025 milliradian error and the lOjt center deter- 
mination error. 

The strongest information source for determining the spacecraft orbital 
parameters is the DSN Doppler data. The actual data are counted Doppler 
frequency over a variable count time. The standard assumption is that the 
Doppler noise is white, which leads to an assumed Doppler error variance 
inversely proportional to the count time. For this analysis, the Doppler 
noise is assumed to be a conservative 1 mm/sec for a one minute count time. 

Filter Algorichm . A sequential filter algorithm was used for this error 
analysis to permit a study of the evolution of the covariance during the close 
encounters. Keplerian orbital elements were taken as the state variables 
for both the spacecraft aM the natural satellite. Tne mass of the particular 
natural satellite under consideration was appended as the thirteenth state 
variable. 

The observational noise discussed above is seldom the limiting error source; 
for satellites of planetary bodies the limiting errors are due to uncertainties 
in the dynamic model. For a Viking spacecraft the dominant errors are t’ie 
higher harmonics of the gravity field and the quasi-stochastic gas leaks from 

7 

the spacecraft . There are a variety of means of including such forces in 

8 9 

the model . For this study, the original Kalman approach was adopted 

because of its simplicity; as will be shown later, this approach is consistent 

with more sophisticated models. 
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The pertinent equations are the observation equations 


Y k “ x k 4 V 


the linearized dynamic model 


*k+l " *k+l, k X k + C k, k+1 W k' 


the mapped covariance 


(-) (♦) T T 

p' ' *. P' ' A 4 H 0 C 

k+1 k+ 1, k k k+1, k k,k+l * k, k+1’ 


the covariance update at the kth observation 


t k* ) * (I - S. V P J'’. 


and the Kalman gain 


•Sc ■ •*'’ \ T > H k p i*’ «*’ ♦ V 1 ' 


The dynamic noise on the twelve orbital elements wan calculutcd by assuming a 

constant, random acceleration between observations. The standard 

deviation of was varied throughout the orbit and whs calculated to be 

equivalent to a 1JJ error in the acceleration due to planetary oblateness at 

the satellite radius from Mars. Although current gravity f ields 1 * ' °* 1 1 have 

formal statistics somewhat smaller than 1JJ, these fields are based on data 

from Mariner 9. which was somewhat insensitive to the gravity field of the 

northern hemisphere. Since the Viking orbiter has a northern hemisphere 

periapsis location, the 1 % uncertainty was chosen to remain conservative. 

-11 ? 

The state noise acceleration ir. bounded below by 10 km/sec to account 
for quasi-stochastic gas leaks. 

Spacecraft Orbit betf-rmlnatiori Accuracy . The Viking navigation team has 
performed detailed studies of the orbit determination accuracies during the 
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nominal Viking mission period . In order to test the validity of the model 
used for this extended mission study, the navigation team results were 
simulated. Figure 3 shows the re ults of thiB comparison. The steady state 
variation of the standard deviation in the spacecraft orbital elements through- 
out the Viking orbit on July U, 1976, is compared with the uncertainties in the 
orbital elements as predicted by the navigation team. The team results (0) 
are based on batch processing of Doppler data from one hour after periapsiB to 
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Hours Past Periapsia 7A/76 
Figure 3. Comparison of Orbit Determination Results 
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one hour before the next periapsis. From this comparison, it was concluded 
that Kalman filter processing is sufficiently consistent to permit application 
of this filter model to the mass determination problem. The major source of 
positional error is an uncertainty of approximately 0.01 degrees In the nodal 
location in the plane of the sky coordinates. The resultant position error 
is therefore proportional to the distance of the spacecraft from the Mars- 
Earth line. Other positional error components are less than 1 km. 


Natural Satellite Orbit Improvement . The uncertainties in the Keplerian 
elements, given in Table 1, indicate that the mujor positional error in 1976 
will be an along-track error of about 30 km for Phobos and 80 km for Deimos. 
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To reduce this error, it is essential that optical imaging of the natural 
satellites be performed prior to the encounter attempt. Optical tracking of 
Deimos will actually be performed during the approach to Mara. It is expected 
that optical tracking of both satellites will be : .*rformed throughout the 
orbital phase of the nominal mission. 

The sequential filter was used to obtain an estimate of the magnitude of the 
downtrack error after orbital phase imaging. Three imuges were simulated for 
each natural satellite during a single orbit in both July and December 1976. 
The resulting improvements in the along-track error are shown in Table 2. 


Table 2 

EFFECT OF OPBI TAL PHASE OPTICAL DATA ON LONGITUDE ERROR (km) 



No Optical Data 

July 

December 

Phobos 

30 

. 

3.0 

Deimos 

80 

1.8 

1.1 


No significant improvement in the other natural satellite orbital elements 
resulted from the optical data. It is assumed that July imaging has been 
performed. Thus, the subsequent analysis utilizes the satellite downtrack 
errors as given by the improved values. 

12 

Analytic Approximation to Mass Determinat ion Accuracy . Anderson introduced 
analytic approximations for the accuracy of determining the mass of an usteroid 
from a single spacecraft flyby. The significant assumptions were (l) that 
the spacecraft trajectory with respect to the asteroid is a rectilinear hyper- 
bola and, (2) that the only trajectory parameter included in the estimation 
process is the distance of closest approach. Eased on these assumptions the 
variance of the estimate of the mass parameter is 

O 2 (GM sin i) « l6bv 3 TO 2 (p)/ir (l) 

If it is further assumed that the distance of closest approach is known, 
the estimate variance becomes 

a 2 (GM sin 1) * 2bv 3 TO 2 (p)/r 


(?) 








10 


In the above formulas, GM is the gravitational constant of the asteroid, i 
is the inclination of the flyby hyperbola to the plane of the sky, b is the 
distance of closest approach, v is the relative velocity, T is the Doppler 
count time, o (p) is the variance on the Doppler velocity measurement. Under 
the basic assumptions of Anderson's analysis, it is impossible to separate 
GM and sin 1, the well-known spectroscopic binary star result. 

The Viking encounters with PhoboB and Deimoo occur over Buch u short time span 
that the above approximations are applicable, producing insight into the 
optimal geometry and expected results. It should be noted that the relative 
velocity is essentially invariant for the degrees of freedom availublt in 
designing the encounter sequence*. Also, the assumption on the Doppler noise 
makes the term W*<p) a constant. These parameters are therefore not available 
for optimizing the experiments. However, assuming thut the inclination to 
the plane of the sky can be determined from the known orbits of the satellites, 
it 1 8 clear that GM is best determined from close encounters with 90° 
inclinations to the plane of the sky. 

Phobos Mass Determination During the January Encounter Owortuni’y . It has 
been shown* that small variations in the phasing and synchronization maneuvers 
permit substantial variations in the encounter geometry. Thus, these maneuvers 
can be used to optimize the encounter sequence to provide the best muss 
determination accuracy. Figure h shows 2 sets of encounters for the January 
Phobos opportunity. These sequences are based on the 1 to 3 synchronized 
orbit3 discussed earlier. Each circle in the figure represents a point where 
the spacecraft descends through the Phobos orbital plane. The components of 
the spacecraft velocity relative to Phobos ure -1.68, 0.^2, and -1.1 km/sec. 

The Earth is about 6° south of the orbital plane in the direction shown. 

Th ! set of encounters represented by open circles (set I) is designed so that 
rich encourter passes nearly through the Phobos-Earth line, thus maximizing 
the inclination to the plane of the sky. As seen in Fig. 5, this design 
produces only two encounters which pass within 100 km of Phobos. The encounter 
set represented by closed circles (set II) is designed to maximize the number 
of close approaches without regard for the inclination to the plane of the 
sky. The intersection points are designed to form a line parallel to the 



Distance of Closest Approach 
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VO Relative 



Figure 1*. Phobos Kncounter Geometry in January 



t . January, 1977 

January, 1977 


Figure 3. Ph>bo3 Kncounter Parameters in January 
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relative velocity vector. This results In 9 encounters which pubs within 
loO km oi' Photo j. 

Both of these encounter sets were analyzed using the sequential filter 
discussed above. Results of the error analyses for both encounter sets are 
presented in Fig. C. The satellite mass a priori was set equal to the nomir.ul 
estimate in Table 1. Encounter set II clearly provid the quickest improve- 
ment in the knowledge of the satellite mass due to the large number of close 
encounters. However, the final uncertainty is about the some for both nets. 



Januory, 1977 

Figure 6. Multi-Oroit Mass Determination 

Neither set pro /ides any substantial improvement in the Phobos rbitui 
elements. The stairstep nature of the muss improvement indicates that very 
little information on the mass remains in the Doppler data after the space- 
craft has receded from the natural satellite. Simple calculations show that 
for a spacecraft-Phobos distance of over a few thousand kilometers, there is 
no observable direct Doppler shift. Therefore, improvements in the mass after 
the encounters must arise from correlations with the spacecraft state variables, 
developed during the close encounter phase. A xack of such improvement could 
be due to the corrupting influence of the Mars gravity field. 
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To confirm this assort ion , the orbit in set II which cane vith'.n 10 km of 
Phobos was analyzed. During the encounter, large correlations are developed 
between the Phobos macs and all six spacecraft elements. After the encounter, 
the general variation of each correlation is fairly well illustrated by 
Fig. 7, in which the correlation between mass and 3<ir.imaJor axis is presented. 
The solid :urve corresponds to the nominal case where the dynamic noise is. 

_j ] 2 

based on a error in J, ar.d a non-gravi tat! oral acceleration of 1C kir./uec . 
The short dashed line corresponds to the case in which the gravity fie’u is 
known perfectly. The third case results from a perfectly known gravity field 



Hours Past Closest Approach 
Figure 7. Evolution of Correlations 


-It P 

and a non-gravitational acceleration of 10 ' km/seo . The uncertainty in the 

gravity field causes the filter to "forget" the encounter well before the next 
periapse passage. Non-gravitational forces nt the 10 ^ levci reduce the 
correlations to .06 at the next ptriapsis. A nearly perfect dynamic model is 
required before information is propagated around to the next encounter. 
Therefore, processing multi-orbits in a continuous mode may provide no 
pivantage over processing single orbits in a multi-batch mode. 


A multi- >atch process was simulated by applying the sequential filter to each 
encounter within set II as if none of the other encounters in that set existed. 
The results of thiB analysis are presented in Fig. 8. Each point represents 
the steady state VLlue of the mass standard deviation ufter the spacecraft 
has receded from Phoboo. The analytical estimates on the standard deviation 
given by equations (i) and (2) above can be directly applied to this filter 
strategy. The Kalman results arc seen to be between the analytic bounds, 
and the general shape of the curves agree well. Thus, the analytic formulas 
provide a computationally efficient tool for future optimization of the 
encounter sequence. 



January, 15.77 

>'igure 8. Gingl e-Orbit Mann Determination 

Statistically combining these individual KnLman covariances in u sequential 
manner provides the lower curve in Fig. 6. A comparison of this result with 
Fig. 6 confirms that a multibatch process in essentially as accurate as the 
single long arc, multi-orbit process. Gince the multi-batch strategy permits 
further optimization with respect to filter parameters, this approach will 
probably be the optimal method for the January encounter. 
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Phobos March Encounter . The March encounter sequence with Phobos vub designed 
to maximize the number of close encounters. The locations of the intersection 
points with the Phobos orbital plane are shown in Fig. 9» there are 9 
encounters within 100 Km of Phobos. Figure 10 shows u large variation in the 
inclination to the plane of the sky. 

The evolution of the mass standard deviation is presented in Fig. 11. The 
final standard deviate 1. January is 3/? larger than the March result. This 
slight improvement is brought nbout by the chang? in the location of the Earth. 
Recall that the major uncertainty in the spacecraft posit ion is produced by a 
rotation of the orbit about the Mars-Earth line. During the January encounters 
the spacecraft is about 9,000 km from this line, whereas in March this distance 
is only U,U0G km. Thus, the relative position of Phobos and the spacecraft in 
slightly better determined in March, which permits the improvement. 


March 23 


VO Relative 
Velocity 


Pho'co 3 Orbital 


Motion 


■*-0 


100 km 


Earth 


• March l6 


Mars 


Figure 9- Phobos Encounter Geometry in March 

Deimos Encounter . It is assumed that the spacecraft has been maneuvered into 
a 1-to-l synchronization with Deimos to provide the maximum number of close 
encounters. This is essential, since the Deimos encounters occur at over 
three times the distance from Mars as the Phobos encounters. This causts 
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March, 1977 

Figure 10. Phobos Encounter Parameters in March 



March, 1977 April 

Figure.* 11. Multi-Orbit Mass Determination in March 
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Plane of Sky, 
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the spacecraft Intersection points within the Deimos orbital plant to neve 
over three tines farther radially per spacecraft orbit. The number of close 
encounters will therefore be limited, even in the optimal case which is 
ilJ"strated in Fig. 12. There are only 5 encounters within 100 km as 
illustrated in Fig. 13. The closest encounter, at 30 km, has a relatively 
high inclination to the plane of the sky. In addition, thic encounter is 
only 7500 km off the Mars-Eart.h line, even though the Deimos encounter phase 
occurs at 23,000 km from Mars. Thus, the relative positions of the spacecraft 
and Deimos are better deter'M. ed than those of the January-Phobos encounter. 


The evolution of the Deimos muss standard deviation is presented in Fig. lU. 
The final uncertainty is 25/f of the a priori estimate of th*. total mnss. 
Almost till of the improvement comes from the encounters on Nov. 29 and 30. 

To determine the contribution of each encounter, the closest encounters in 
this sequence were reanalyzed in the single encounter mode. These resuits 
are also presented in Fig. ]b. These are the only two encounters which 
provide better than $0% mass determinat ion. 


Deimos 

Orbital 


Motion 


Nov. 31 • 




T Earth 


200 km 



Mars 




Nov. 27 • 


\ 



V0 Relative 
Velocity 


Figure 12. Deimos Encounter* Geometry 


The above analysis is based on the 1-to-l synchronization between orbital 
periods, vhich may require a propulsive expenditure of about 50 m/s. A 
h-to-5 orbital period ratio can be established for less than 5 m/s. However, 
such suboptimal encounter sequences provide only one close encounter, limiting 
the accuracy of the mass determination to about 35$. In addition, these 
sequences do not provide independent estimates of the mass. 
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Figure 13. Deimos Encounter Parameters 
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Probability of Collision v 1th Phobos and Dtlnos . Determining the mass of the 
natural satellite depends strongly on the minimum distance of closest approach. 
This distance must be sufficiently large to assure a low probability of 
spacecraft impac* . The impact probability has been determined by assuming 
that the encounte. • ometries are those diocucted earlier. Based on the 
spacecraft u*d satt ite orbit determination accuracies at closest approach 
minus 1 , 2 , and 5 orbits, the probability of impact during the closest 
encounter was calculated using n Monte Carlo technique with 10,000 samples. 

The results are presented in Table 3. 


Table 3 


Percent Impacts 

Orbits 

to Closest Approach 

5 

2 

1 

J anuary 

7.? 

7.8 

1 . 1 * 

March 

?.U 

0.7 

0.1 

November 

3.1 

X.k 

1.0 


To assure an impact probability of less than 1 # will require a decision one 
orbit prior to the closest encounter except in March, when a reduced collision 
probability results from the previously discussed improvement, in spacecraft 
determination. 

A simulation was also performed for the January case excluding the downtrack 
improvement afforded by the optical data. The probability of collision was 
increased from 1.1»J{ to 3.8?. Thus, optical data plays an important role in 
the design of close encounters. 

CONCLUDING REMARKS 

The application of a Kalman filter to the l'h^bos and Deiraos encounter 
oppor. unities during the Viking extended mission has shown that the mass of 
Phobos can be determined to the 10? level, during either the January or 
March (1977) opportunities. The Deimos mass can be determined to about L’5? 
with the optimal encounter geometry, which provides two close encounters. 

This uncertainty is Increased to about. 35?, if only one close encounter is 
possible. 
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T 
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M 


Si ♦ 1, k 

R 
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v 

W k 

Si 
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U 


aemimajor axis of Keplerion orbit, kin 
distance of closest approach, km 
state noise transition matrix 

eccentricity of Keple-ian orbit 
observation noise vector 

product of gravitational constant and mass of central boay, 
km^/sec^ 

altitude of apoapsis, km 
altitude of periapsis, km 

observation matrix 

orbital inclination, deg 
identity matrix 

planetary oblatencss coefficient 

Kalman gain matrix 

mean anomaly 

mean motion, sec J 
covariance matrix 

state noise covariance matrix 

mean radius of Mars, km 
observational noise covariance matrix 

relative velocity, km/sec 
random acceleration 

state vector 

observational vector 

variance on parameter x 

2 

variance on the Doppler velocity measurement, (kro/sec) 
Doppler count time, sec 
argument of periapsis, deg 
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? 1 


n longitude of ascending node of Keplerlun orbit, deg 

Subscripts: 

k kth observation 

Superscripts: 

T transpose 

-1 Inverse 

(-) prior estimate 

(♦) updated estimate 

A dot over symbol denotes differentiation with respect to time. 
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